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Using density functional theory, we study physical properties of boron nitride nanotubes (BNNTs) 
with the substitutional carbon pair defect. We also consider the Stone- Wales (SW) rearrangement 
of the C-C pair defect in the BNNT. The formation energy of an SW defect of the carbon dimer 
is approximately 3.1 eV lower than that of the SW-transformed B-N pair in the undoped BNNT. 
The activation energies show that the SW defect in the C-doped BNNT may be experimentally 
observed with a higher probability than in the undoped BNNT. Finally, we discuss the localized 
states originating from the carbon pair impurities. 



As a wide-band-gap nanomaterial[ll-Q, hexagonal 
boron nitride (h-BN) is a stable crystalline form con- 
sisting of equal numbers of boron and nitrogen atoms in 
a honeycomb arrangement UQ- As it has considerably 
high thermal and chemical stabilities 0, it can be used 
for high-temperature applications. The BN nanotube 
(BNNT) [1, Q is a nano-sized seamless cylinder that can 
be considered a rollcd-up h-BN sheet. It also has a wide 
band gap that is almost independent of the tube diame- 
ter, and chiralitv[lo|. For sp^ -bonded carbon nanostruc- 
tures such as fullerenes, carbon nanotubes and graphene, 
the Stone- Wales (SW) transformationfllj is believed to 
introduce topological defects or isomerization, which re- 
sults in a 90° rotation of two carbon atoms with respect 
to the midpoint of the C-C bond. 

When an SW defect is present in h-BN or BNNT, B- 
B and N-N bonds should be created, and the formation 
of these bonds increases the total energy of the system. 
Thus, the topological defect created by the SW rear- 
rangement is lacking. However, the formation of B-B 
and N-N bond pairs can be avoided if there is a C-C pair 
defect in h-BN (or BNNT). Using post-synthesis dop- 
ing, carbon impurities can be substituted in the h-BN 
lavers jl2Ml4l |. During the formation of BNNTs and h- 
BN sheets using the chemical vapor deposition (CVD), 
carbon atoms may act as substitutional defects, because 
residual hydrocarbon may remain in the chamber, or car- 
bon impurities are solubilized in the metal catalyst. If 
carbon substitutional atoms congregate together in a lo- 
cal area in h-BN (or a bundle of BNNTs) by segrega- 
tion, they may form the interlayer (or intertube) con- 
duction channels. As mentioned above, the C-C pair de- 
fect in the h-BN network prohibits the formation of B-B 
and N-N bonds during the SW transformation. As the 
SW-transformed defects may affect the mechanical and 
electrical properties of the BN nanostrucutres, we should 
thus study these defects in the C-doped BN sheets or BN- 
NTs. This work may explain the formation of topological 
defects in the carbon-doped BNNTs or h-BN layers that 
can be produced in the CVD processes. In this paper, 
we present our first-principles study of the structural and 



electronic properties of the BNNTs containing two sub- 
stitutional carbon defects. The effect of carbon doping 
on the SW transformation is also examined. 

We carried out total energy calculations for our model 
systems based on the density functional theory 15| . The 
local density approximation (LDA)[l6| with spin polar- 
ization was used for the exchange-correlation functional. 
Some calculations were repeated using the generalized 
gradient approximation. We found that the exchange- 
correlation functional did not change our main conclusion 
that the SW defect in the C-doped BNNT may be ob- 
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FIG. 1: Optimized structures of (8,0) BNNTs with a C-C 
pair substitutional defect or the SW defects, (a) The BNNT 
with a C-C pair substitutional defect parallel to the tube axis, 
(b) the BNNT with a C-C pair substitutional defect with 
an oblique angle with respect to the tube axis, (c) and (d) 
show the SW-transformed C-C pair defect from (a) and (b), 
respectively, (e) and (f) are the BNNTs with SW defects 
which are parallel and tilted to the tube axis, respectively. 
Yellow, green and black balls represent nitrogen, boron and 
carbon atoms, respectively. 
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served with a higher probabihty than in the bare BNNT. 
To study the energetics and the activation energy bar- 
rier in the SW transformation of the systems, we used a 
plane wave basis set with an energy cutoff of 400 eV 
and ultrasoft pseudopotentials[l3l implemented in the 
Vienna ab initio simulation package (VASP)[l3,[lll- The 
atomistic models were relaxed until the residual forces on 
the atoms became smaller than 0.03 eV/A. For the ac- 
tivation ener gy b arrier, we employed the nudged elastic 
band mcthod[20j. Eight replicas were selected, including 
the initial and final configurations to construct an clastic 
band. To plot the density of states (DOS) and the wave- 
functions with phase information, we used the OpcnMX 
code2l|. In the calculations using the OpenMX code. 



the ionic potentials were described with norm-conserving 
TrouUier-Martins pseudopotentials. Wavefunctions were 
expanded in a pseudoatomic orbital basis set 22[ with an 
energy cutoff of 120 Ry. As a model system, we chose 
the (8,0) zigzag BNNT containing 48 B atoms and 48 N 
atoms. The supercell size in the lateral direction was 25 
A to avoid the interaction between the defects in neigh- 
boring BNNTs, and the size in the axial direction was 
12.96 A. 




Reaction coordinate 




First, we optimized the C-doped and undoped BN- 
NTs. The perfect BNNT has a B-N bond length of - 
1.44 A. Our total energy calculations show that the two 
substitutional carbon atoms tend to bind to each other. 
This finding agrees with K. Yuge's report that the lowest- 
formation-energy structures for monolayer boron-carbon- 
nitridc favor B-N and C-C bonds whereas they disfavor 
B-C, C-N, B-B, and N-N bondsjH]. For C-doping, two 
types of C-C pairs can be considered. As shown in Figure 
1, one pair (BNCl) was parallel to the tube axis and the 
other (BNC2) had an oblique angle with the tube axis. 
The bond length of the C-C pair was 1.37 A (1.39 A) for 
BNCl (BNC2). Note that the bond length of benzene 
(CgHg) has a very similar C-C bond length of 1.40 A. Al- 
though the model of BNCl had the lowest energy among 
our model systems for the C-C pair, the energy differ- 
ence between BNCl and BNC2 was only ~1 meV/atom. 
Thus, two structures may coexist. In fact, by the chem- 
ical potential analysis, we found that the formation of 
carbon substitutional defects in the BNNTs, BNCl or 
BNC2, is possible when the chemical potential of an iso- 
lated C-C pair is larger than by -18.02 cV or -17.90 eV. 
Note that the energy of an C-C pair is -10.09 eV. 

In the C-C pair, the C atom bonded to two B atoms 
protruded out by ^ 0.1 A. The B-C bond length was 1.52 
A and the N-C bond length was 1.46 A, which was due to 
the longer covalent radius of B (jcovaient,B = 0.84 A) than 
that of N (rcot.a/ent.w = 0.71 A). When the C-C bond was 
rotated by 90°, the bond length slightly decreased to 1.35 
A. Here, the B-C bond length was increased to 1.61 A, 
and the N-C bond length was slightly decreased to 1.44 
A. After the SW transformation from the BNCl to the 
BNCRl structure, the total energy of BNCRl increased 
by 2.5 cV, compared to that of BNCl. For this structure, 
the activation energy from BNCl to BNCRl was 5.5 eV. 
Similarly, the formation and activation energies were 2.2 
eV and 5.6 eV, respectively, for the SW transformation 
from BNC2 to BNCR2 (see Table I). 

For comparison, we also calculated the formation and 
activation energies for the SW transformation in the un- 
doped (8,0) BNNT as shown in Figure 2. For the SW 
transformation from BNl to BNRl (or from BN2 to 
BNR2), the formation energy was 5.6 eV (5.7 cV), and 




FIG. 2: (Top panel) Energy profile along the reaction co- 
ordinate in the SW transformation of BNNTs with a C-C 
pair substitutional defect or the SW defects. The states "A", 
"B", and "C" represent the initial, transition, and final states 
in the SW transformations, respectively. The activation and 
formation energies are denoted as Ea and Er, respectively. 
(Bottom panel) The corresponding three states are shown for 
the SW transformation from BNCl to BNCRl. 



TABLE I: Calculated activation {Ea) and formation {Er) en- 
ergies for the SW transformation from the C-doped or the 
BNNTs to the BNNTs with SW defects. 



Ea (eV) Er (eV) 



BNCl - 


■> BNCRl 


5.50 


2.51 


BNC2 - 


■> BNCR2 


5.64 


2.17 


BNl - 


^ BNRl 


7.04 


5.61 


BN2 - 


^ BNR2 


7.63 


5.72 
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the activation energy was 7.0 eV (7.6 eV). Here, a B-B 
bond and a N-N bond are produced, and their lengths are 
1.73 A (B-B) and 1.47 A (N-N), respectively. We found 
that the rotating C-C bond (1.24 A) was shorter than the 
rotating B-N bond (1.29 A) in the transition states. This 
implies that the shorter bond results in a lower activation 
energy barrier during the SW transformation in a small 
area. The high formation energy is due to the unstable 
B-B and N-N bonds as mentioned earlier. Interestingly, 
electrons were accumulated in the N atom of the rotating 
B-N pair in the transition state for the undoped BNNT. 
In addition, our charge analysis revealed that electron do- 
nation from the C-C pair defect in the C-dopcd BNNT 
slightly weakens the dangling-bond-like character of the 
B and N atoms in the transition configuration. In our 
previous study of the SW transformation in M@C6o (M 
= K, Ca, and La), we similarly observed that the elec- 
tron donation by the incorporated metal atom lowers the 
activation energy barrier [24|. These findings can explain 
the lower activation energy barrier of the C-C pair defect 
compared to that of the B-N pair. Therefore, when the 
C-dopcd BNNT does exist, e.g., from the sample prepa- 
ration, we can conclude that the probability that the SW 
defect is observed is higher in the C-doped BNNT than 
in the undoped BNNT. 

To check a charge doping effect, we calculated doping 
of a hole (charge +e state) for all the model structures 
(BNRl, BNR2, BNCRl, and BNCR2), and an electron 
doping case (charge -e state) for BNR2 structure. When 
a hole was doped, the activation energy barriers were 
reduced from 7.04 eV to 6.91 eV for BNRl, 7.63 eV to 7.18 
eV for BNR2, 5.50 eV to 5.23 eV for BNCRl, and 5.64 
eV to 5.52 eV for BNCR2, respectively. For an electron 
doping in BNR2, and the activation energy was decreased 
from 7.63 eV to 6.69 eV. Therefore, we can conclude that 
charge doping affects the activation energy barriers of the 
SW transformation. 

In a growth process of the BNNTs, carbon atoms may 
be incorporated as residual impurities. Or one can pro- 
duce carbon substitutional defects in the BNNTs on pur- 
pose for doping. Our results demonstrate that two car- 
bon atoms far from each other tend to become closer and 
form a C-C pair. The energy difference between the C-C 
pair and two separate carbon atoms was larger than 5 eV 
in our LDA calculation. On the other hand, some metal 
atoms or molecules such as O2 and Br2 may be adsorbed 
in the vicinity of the C-C pair. Then, electron transfer 
takes place between the adsorbate and the nanotube. As 
mentioned above, charge doping enables the activation 
energy barrier to become lowered and the probability of 
the SW transformation to become higher. 

According to our LDA calculation, the energy band 
gap of the perfect (8,0) BNNT is 3.6 eV as shown in Fig- 
ure 3a. R. Wu et al. reported that the BNNT could 
have a magnetic moment due to the carbon doping [25|. 
They considered only one carbon atom in the superccll. 
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FIG. 3: (a) Total DOS of the perfect (8,0) BNNT. Total DOS 
of BNNTs with (b) the C-C pair substitutional defect or (d) 
the SW defects. Dotted lines in (c) and (e) indicate the pro- 
jected DOS of two C atoms in the C-C pair defect in each 
case. 



However, our calculational results demonstrated that the 
magnetic moment does not occur in the BNNT with the 
C-C pair defect. According to a recent report by N. 
Berseneva et al., the magnetic moment is zero when the 
numbers of carbon atoms replaced by nitrogen and boron 
atoms are identical (iij. Figure [3}d shows the total DOS 
of BNNTs with the C-C pair substitutional defect (e.g., 
BNCl), which clearly demonstrates that a few localized 
states occur in the band gap in the presence of the C-C 
pair defect. As shown in Figure the projected DOS 
of the two C atoms in the C-C pair defect indicates that 
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the gap states originate from the substitutional carbon 
pair defects. When a C-C pair forms a SW defect by 
a 90° rotation as seen in BNCRl, the electronic states 
originating from the SW defect have flatter bands with 
more localized characters. These characteristics are re- 
flected in the sharp and narrow peaks in the DOS in the 
gap (Figure [3jl and[3j5). 

Figure 2] shows wavcfunction characters resulting from 
the defect states in the BNNTs with the substitutional 
carbon pair defect and its SW-transformcd defect. The 
occupied state at —1.4 cV and the unoccupied state at 
+1.7 cV shown in Figure|4^ are the localized defect states 
originating primarily from the C-C pair of the BNCl 
structure. The filled state has a bonding character, 
whereas the unfilled state has a antibonding character 
at the C-C pair. As displayed in Figure |3^, the electron 
density at the C atom bonded to B atoms is somewhat 
higher than that at the C atom bonded to N atoms at 
— 1.4 eV. Additionally, the B-C bonds are in the bonding 
states but the C-N bonds are in the antibonding states. 
On the other hand. Figure IUd shows the localized defect 
states at —1.3 cV and +1.8 cV from the SW-transformed 
C-C pair defect of the BNCRl structure. Similar to the 
BNCl case, bonding and antibonding characters of the 
carbon pair appear at —1.3 eV and +1.8 eV, respectively. 
Note that there is no mirror (reflection) symmentry for 
the BNNT with respect to the plane containing the mid- 
point of the C-C bond. However, the mirror symmetry 
is manifest with respect to the planes containing the C- 
C bond itself (as in BNCl) as well as the midpoint of 
the SW-transformed C-C bond (as in BNCRl). These 
features were retained for two localized states as shown 
in Figure 11^ and also for a localized state at —1.3 eV as 
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FIG. 4: Wavefunctions originating from the defect states 
in the BNNTs with (a) the substitutional C-C pair de- 
fect (BNCl) and (b) the SW-transformed C-C pair defect 
(BNCRl). Red and blue colors represent the opposite signs 
in the wavefunctions at the F-point. 



shown in Figure HJd. 

In summary, we have introduced substitutional carbon 
pair defects in BNNTs to examine the modifications in 
the structural and electronic properties of the nanotube. 
Our results demonstrate that the two carbon atoms pre- 
fer to bind to each other with a bond length of ^1.37 
A. For the SW transformation of the C-C pair defect, 
the formation energy of the SW defect is approximately 
2.5 eV and the activation energy is approximately 5.5 
eV. Both the formation and activation energies of the 
C-C bond are much lower than those of the B-N bond 
in the undoped BNNT in the SW transformation, which 
are ~5.5 eV and ~7.0 eV, respectively. Consequently, it 
would be very difficult to find a topological defect such as 
the SW defect in the perfect BNNT. Conversely, such a 
topological defect could be observed with a higher prob- 
ability in the C-rich defective BNNTs. 
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